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A modification to the Priedmann Robertson Walker equation is proposed 
^ . in which the universe is flat, matter dominated, and accelerating. An ad- 
0^ \ ditional term, which contains only matter or radiation (no vacuum cantri- 
ps \ bution), becomes the dominant driver of expansion at a late epoch of the 
T-H ■ universe. During the epoch when the new term dominates, the universe ac- 
' celerates; we call this period of acceleration the Cardassian era. The universe 
. can be flat and yet consist of only matter and radiation, and still be compat- 
\ ible with observations. The energy density required to close the universe is 
much smaller than in a standard cosmology, so that matter can be sufficient 
O ' to provide a flat geometry. The new term required may arise, e.g., as a conse- 
' quence of our observable universe living as a 3-dimensional brane in a higher 
■ dimensional universe. The Cardassian model survives several observational 
^ ! tests, including the cosmic background radiation, the age of the universe, the 
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Recent observations of Type lA Supernovae as well as concordance with other 
observations (including the microwave background and galaxy power spectra) indicate that 
the universe is accelerating. Many authors have explored a cosmological constant, a decaying 
vacuum energy and quintessence [^]-|^ as possible explanations for such an acceleration. 

Here we propose an alternative which invokes no vacuum energy whatsoever. In our 
model the universe is flat and yet consists only of matter and radiation. Pure matter (or 
radiation) alone can drive an accelerated expansion if the Friedmann Robertson Walker 
(FRW) equation is modified by the addition of a new term on the right hand side as follows: 

H^ = Ap + Bp'', (1) 

where H = R/R is the Hubble constant (as a function of time), R is the scale factor of the 
universe, the energy density p contains only ordinary matter and radiation, and we will take 

n < 2/3. (2) 

In the usual FRW equation B = 0. To be consistent with the usual FRW result, we take 
A = We note here that the geometry is fiat, as required by measurements of the cosmic 

background radiation ||^, so that there are no curvature terms in the equation. There is no 
vacuum term in the equation. This paper does not address the cosmological constant (A) 
problem; we simply set A = 0. 

In this paper, we first study the phenomenology of the ansatz in Eq.(|l|), and then turn to 
a discussion of the origin of this equation[|. Directions for a future search for a fundamental 
theory will be discussed. 

The Role of the Cardassian Terrr^: The new term in the equation (the second term on 
the right hand side) is initially negligible. It only comes to dominate recently, at the redshift 
Zeg ~ 0(1) indicated by the supernovae observations. Once the second term dominates, it 
causes the universe to accelerate. We can consider the contribution of ordinary matter, with 

p = PoiR/Ro)-^ (3) 

to this new term. Once the new term dominates the right hand side of the equation, we 
have accelerated expansion. When the new term is so large that the ordinary first term can 
be neglected, we find 

Rocm (4) 



^As discussed below, we were motivated to study an equation of this form by work of Chung and 



Freese [10| who showed that terms of the form can generically appear in the FRW equation as 



a consequence of embedding our observable universe as a brane in extra dimensions. 

^The name Cardassian refers to a humanoid race in Star Trek whose goal is to take over the 
universe, i.e., accelerated expansion. This race looks foreign to us and yet is made entirely of 
matter. 
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so that the expansion is superluminal (accelerated) for n < 2/3. As examples, for n = 2/3 
we have R ^ t; for n = 1/3 we have R ~ t^; and for n = 1/6 we have R ~ t^. The case 
of n = 2/3 produces a term in the FRW equation oc R^^; such a term looks similar to 
a curvature term but is generated here by matter in a universe with a flat geometry. Note 
that for = 1/3 the acceleration is constant, for n > 1/3 the acceleration is diminishing in 
time, while for n < 1/3 the acceleration is increasing (the cosmic jerk). 

The second term starts to dominate at a redshift z^q when Ap{zeq) = Bp'^{zeq), i.e., when 

B/A = pl--{l + Zeqf'-^\ (5) 

From evaluating Eq.(|l]) today, we have 

H', = Apo + Bp-, (6) 

so that 

A = H'Jpo-Bpr'. (7) 
From Eqs.(^ and (0), we have 

pS[l + (l + Ze,)3(l-")]' ^ ^ 

We have two parameters in the model: B and n, or, equivalently, Zeq and n. Note that i? 
here is chosen to make the second term kick in at the right time to explain the observations. 
As yet we have no explanation of the coincidence problem; i.e., we have no explanation for 
the timing of Zf.q. Such an explanation would arise if we had a reason for the required mass 
scale of B. The parameter B has units of mass^~^". Later, we will discuss the origin of the 
Cardassian term in terms of extra dimensions, and discuss the origin of the mass scale of 
B. As discussed below, to match the CMB and supernovae data we take 0.3 < Zeq < 1, but 
this value can easily be refined to better fit upcoming observations. 
What is the Current Energy Density of the Universe? 

Observations of the cosmic background radiation show that the geometry of the universe 
is flat with = 1. In the Cardassian model we need to revisit the question of what value 
of energy density today, Pq, corresponds to a fiat geometry. We will show that the energy 
density required to close the universe is much smaller than in a standard cosmology, so that 
matter can be sufficient to provide a flat geometry. 

The energy density po that satisfles Eq.(|) is, by deflnition, the critical density. From 
Eqs.([|) and (|^), we can write 

H' = A[p + pl-'-{l + z,,f'~-^p-]. (9) 

Evaluating this equation today with A = Sn / {3mpi) , we have 

^o = ^po[^+a+^eq?^'-''^ (10) 
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Defining po = ^oPc we find that the critical density pc has been modified from its usual 
value, i.e., the number has changed. We find 

- - ° ^' (11) 



87r[l + (l + ze,)3(i-)]" 
Thus 

pc = Pcoid X F{n) (12) 

where 

F(n) = [l + (l + ;.,,)3(l-")]-l (13) 

and 

p,,oid = 1.88 X 10-29/i2gni/cm-=^ (14) 

and ho is the Hubble constant today in units of 100 km/s/Mpc. In Figure 1, we have plotted 
the new critical density function of the two parameters n and Zeq. For example, if 

we take Zeq = 1, we find 

F = (1/3, 1/5, 0.15) for n = (2/3, 1/3, 1/6) respectively. (15) 

We see that the value of the critical density can be much lower than previously estimated. 
Since = 1, we have today's energy density as po = pc as given above 0, i.e., 

Po = (1/3, 1/5,0.15) X 1.88 X 10"^^/t^ gm/cm^ for n = (2/3, 1/3, 1/6) and = 1- (16) 

For larger values of Zeq, the modification to the value of pc can be even larger. Note the 
amusing result that for Zeq = 2 and n = 1/12, we have pc = 0.046pc,oZd so that baryons would 
close the universe (not a universe we advocate). 

Cluster Baryon Fraction 

For the past ten years, a multitude of observations has pointed towards a value of the 
matter density Po ~ Q-^Pcoid- The cluster baryon fraction [|11],|12[ as well as the observed 
galaxy power spectrum are best fit if the matter density is 0.3 of the old critical density. 



^An alternate possible definition would be to keep the standard value of pc and discuss the 
contribution to it from the two terms on the right hand side of the modified FRW equation. Then 
there would be a contribution to 17 from the p term and another contribution from the term, 
with the two terms adding to 1. This is the approach taken when one discusses a cosmological 
constant in lieu of our second term. However, the situation here is different in that we have only 
matter in the equation. The disadvantage of this second choice of definitions would be that a value 
of the energy density today equal to pc according to this second definition would not correspond 
to a flat geometry. 
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Recent results from the CMB also obtain this value. In the standard cosmology this 
result implied that matter could not provide the entire closure density. Here, on the other 
hand, the value of the critical density can be much lower than previously estimated. Hence 
the cluster motivated value for is now compatible with a closure density of matter, ^0 = 1, 
all in the form of matter. For example, if n = 0.6 with Zeq = 1, or if n = 0.2 with Zeq = 0.4, 
a critical density of matter corresponds to po ~ 0.3pc,o«d, as required by the cluster baryon 
fraction and other data. In Figure 1, one can see which combination of values of n and Zgq 
produce the required factor of 0.3. If we assume that the value po = 0.3pc,o«d is correct, for 
a given value of n (that is constant in time) we can compute the value of z^q for our model 
from Eq. (|I3|) . Table I lists these values of n and Zeq- Henceforth, we shall focus on these 
combinations of parameters. 
Age of the Universe 

In the Cardassian model, the universe is older due to the presence of the second term. 
In Table I, we present the age of the universe for various values of n (under the assumption 
that Po = O.Spcoid)- 

If one takes to > lOGyr as the lower bound on globular cluster ages, then one requires 
toHo > 0.66 for ho = 0.65. If one requires globular cluster ages greater than 11 Gyr p3| , 
then toHQ > 0.73 for = 0.65. All values in Table I satisfy these bounds. 

Structure Formation 

Since the new (Cardassian) term becomes important only at z^q ~ 1, early structure 
formation is not affected. Below we discuss the impact on late structure formation during 
the era where the Cardassian term is important. This term accelerates the expansion of the 
universe, and freezes out perturbation growth once it dominates (much like when a curva- 
ture term dominates); this freezeout happens late enough that it is relatively unimportant. 
To obtain an idea of the type of effects that we may find, instead of analyzing the exact 
perturbation equations with metric perturbations included, we will merely modify the time 
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_n Zeq Hptp 

0.60 1.00 0.73 

0.50 0.76 0.78 

0.40 0.60 0.83 

0.30 0.50 0.87 

0.20 0.42 0.92 

0.10 0.37 0.95 

0.00 0.33 0.99 



TABLE I. Values of Zeq for various values of n corresponding to a universe with pQ = O.Spc^oid- 
The age of the universe today to corresponding to the two parameters n and Zgg is listed in the 
last column, where Hq is the value of the Hubble constant today. 

dependence of the scale factor in the usual Jeans analysis equation. For now we take the 
standard equation for perturbation growth; as a caveat, we warn that recent structure for- 
mation may be further modified due to a change in Poisson's equations as described below. 
For we now we take 

6 + 2{R/R)6 = inpS/mli (17) 

where 6 = {p — p)/p is the fluid overdensity. Now one must substitute Eq.(|l|) for R/R. In 
the standard FRW cosmology with matter domination, R ~ t^/^, and there is one growing 
solution to S with 5 ^ R ^ t^^^. This standard result still applies throughout most of the 
(matter dominated) history of the universe in our new model, so that structure forms in the 
usual way. 

Modifications set in once the new Cardassian term becomes important. When R ~ t^, 
Eq.(|17D can be written 

S"{x) + '^S'ix) - ^-x-^n = 0, (18) 

where x = t/tQ with to denoting the time today and superscript prime refers to d/dx. This 
equation can generally be solved in terms of Bessel functions for constant p (such as is the 
case once the Cardassian term completely overrides the old term). A simple example is the 
case of n = 2/3 and p = 1; in the limit x » 3/4, the last term in Eq.([TB|) can be dropped 
and the equation is solved as S{t) = ai + a2t~^. Perturbations cease growing and become 
frozen in. This result agrees with the expectation that in a universe that is expanding more 
rapidly, the overdensity will grow more slowly with the scale factor. As mentioned at the 
outset, as long as the Cardassian term becomes important only very late in the history of the 
universe, much of the structure we see will have already formed and be unaffected. Further 
comments on late structure formation (e.g. cluster abundances) follow below. 
Doppler Peak in Cosmic Background Radiation 

Here we argue that the location of the first Doppler peak is only mildly affected by 
the new Cardassian cosmology. We need to calculate the angle subtended by the sound 
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horizon at recombination. In the standard FRW cosmology with flat geometry, this value 
corresponds to a spherical harmonic with / = 200. A peak at this angular scale has indeed 
been conflrmed p]. In the Cardassian cosmology we still have a flat geometry. Hence, we 
can still write 

e = sjd, (19) 

where s,, is the sound horizon at the time of recombination tr and d is the distance a light 
ray travels from recombination to today. To calculate these lengths, we use the fact that for 
a light ray ds"^ = = —dt"^ + a^dx^ to write 

d= f° dt/a. (20) 

Jtr 

Following the notation of Peebles we deflne the redshift dependence of H as 

H{z) = HoE{z) (21) 



so that Eq.(l^) can be written 



d = / — -. (22) 



Similarly, the sound horizon at recombination is 

POO 

= / dt/a. (23) 



In standard matter dominated FRW cosmology with fim,o = 1, E{z) = (1 + 2;)^/^ in 

Eq.(|2D and rf = 2/i^o^o- 

For the cosmology of Eq. (|l|) , we have 



E{zY = F X (1 + z)^ + (1 - F) X (1 + zY' (24) 

with F given in Eq.(|l^). As discussed previously, as our standard case we will take F = 
Pc/ Pc,oid = 0.3. With this assumption, and by using expression Eq.(piD in Eq. P^) , we flnd 
that d changes by a factor of (1.47, 1.88, 2.04, and 2.23) for n =(0.6, 0.3, 0.2, and 0.1) 
respectively compared to the the usual (nonCardassian) FRW universe with po = pc,oid- In 
addition cx -^In^SKjh^^rill ^here r, = 0.042(F/i2)-i and = SO^bh^ and we use 
h = 0.7 and ^f, = 0.04. We' flnd that changes by a factor of (1.44, 1.62, 1.67, 1.29) for 
n=(0.6, 0.3, 0.2, and 0.1) respectively compared to the usual FRW universe with po = Pc,oid- 
The angle subtended by the sound horizon on the surface of last scattering decreases and 
the location (/) of the flrst Doppler peak increases by roughly a factor of 

(1.02, 1.11, 1.12, 1.13) for n= (0.6, 0.3, 0.2, 0.1) respectively (25) 

compared to the usual FRW universe with po = Pc,oid- This shift still lies within the experi- 
mental uncertainty on measurements of the location of the Doppler peak. 
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We note the following: in the same way that a nonzero A may make the current CBR 
observations compatible with a small but nonzero curvature, indeed a nonzero Cardassian 
term could also allow for a nonzero curvature in the data. A more accurate study of the 
effects of Cardassian expansion on the cosmic background radiation (including the first and 
higher peaks) is the subject of a future study. 

The Cutoff Energy Density 

An alternate way to write Eq.(|ip is 

H' = Ap[l + ip/pcutoffr-'] (26) 

where Pcutoff = p{zeq) = A/ B. This notation offers a new interpretation; it indicates that 
the second term only becomes important once the energy density of the universe drops below 
Pcutoff, which has a value a few times the critical density. Hence, regions of the universe 
where the density exceeds this cutoff density will not experience effects associated with the 
Cardassian term. In particular, we can be reassured that the new term won't affect gravity 
on the Earth or the Solar System. The density of water on the Earth is 1 gm/cm^, which is 
28 orders of magnitude higher than the critical density. 
Comparing to Quintessence 

We note that, with regard to observational tests, one can make a correspondence between 
the Cardassian and Quintessence models; we stress, however, that the two models are entirely 
different. Quintessence requires a dark energy component with a specific equation of state 
[p = wp), whereas the only ingredients in the Cardassian model are ordinary matter {p = 0) 
and radiation [p = 1/3). However, as far as any observation that involves only R{t), or 
equivalently H{z), the two models predict the same effects on the observation. Regarding 
such observations, we can make the following identifications between the Cardassian and 
quintessence models: n ^ w + 1, F ^ Qm, and 1 — F ^ Qq, where w is the quintessence 
equation of state parameter, Qm = Pm/ Pc,oid is the ratio of matter density to the (old) 
critical density in the standard FRW cosmology appropriate to quintessence, Qq = PqI Pc,oid 
is the ratio of quintessence energy density to the (old) critical density, and F is given by 
Eq. (p!3D . In this way, the Cardassian model can make contact with quintessence with regard 
to observational tests. 

All observational constraints on quintessence that depend only on the scale factor, R{t) 
(or, equivalently, H[z)) can also be used to constrain the Cardassian model. However, 
because the Cardassian model requires modified Einstein equations (see below), the gravita- 
tional Poisson's equations and consequently late-time structure formation may be changed; 
e.g., the redshift dependence of cluster abundance should be different in the two models. 
These effects (and others, such as the fact that quintessence clumps) may serve to distin- 
guish the Cardassian and quintessence models. The correspondence with quintessence, as 
well as discussion of distinguishing tests will be the subject of a future paper. 

Best Fit of Parameters to Current Data 

We can find the best fit of the Cardassian parameters n and z^q to current CMB and 
Supernova data. The current best fit is obtained for po = 0.3pc,oid (as we have discussed 
above) and n < 0.4 (equivalently, w < —0.6) ppj , !^ . In Table I one can see the values 
of compatible with this bound, as well as the resultant age of the universe. As an 
example, for n = 0.2 (equivalently, w = —0.8), we find that z^q = 0.42. Then the position 
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of the first Doppler peak is shifted by a factor of 1.12. The age of the universe is 13 Gyr. 
The cutoff energy density is Pcutoff = 2.7pc, so that the new term is important only for 
P < Pcutoff = '^■7pc- Hence, as mentioned above, the Cardassian term won't affect the 
physics of the Earth or solar system in any way. 

We note the enormous uncertainty in the current data; future experiments (such as SNAP 



2^ ) will constrain these parameters further. 
Extra Dimensions 

A Cardassian term may arise as a consequence of embedding our observable universe 
as a 3+1 dimensional brane in extra dimensions. Chung and Freese showed that, in a 
5-dimensional universe with metric 

ds"^ = -q^{T, u)dT'^ + a^(r, u)dx^ + 6^(r, u)du^, (27) 

where u is the coordinate of the fifth dimension, one may obtain a modified FRW equation 
on our observable brane with H"^ ~ p" for any n (see also |19[)- This result was obtained with 



5-dimensional Einstein equations plus the Israel boundary conditions relating the energy- 
momentum on our brane to the derivatives of the metric in the bulk. 

We do not yet have a fundamental higher dimensional theory, i.e., a higher dimensional 
T^u, which we believe describes our universe. Once we have this, we can write down the mod- 
ified four-dimensional Einstein's equations and compute the modified Poisson's equations, 
as would be required, e.g., to fully understand latetime structure formation. 

There is no unique 5-dimensional energy momentum tensor T^y that gives rise to Eq. (|1]) 
on our brane. Hence, in this paper we construct an example which is easy to find but 
is clearly not our universe, simply as a proof that such an example can be written down. 
Following [|lOl (see Eqs. (24) and (25) there with F{u) = m), we have constructed an example 
of a bulk T^y for arbitrary n in ~ p", matter on the brane as in Eq.(^), and with q = h 
in Eq. (|27|) . We display only Tq here (the other components will be published in a future 
paper): 



2rpQ 



3 n B ^ € 



5-^0 ~ 2 2 



where 



exp 



4 • 81"Sn - 16"/tf(— )" (nV^ + 

UT 



, . 2 + n „ 1 r, , \ , 2 

i2/3) — B-^Kl{—)nu 



(28) 



(29) 



and the constant is related to the 5-dimensional Newton's constant G5 and 5-D reduced 
Planck mass M5 by the relation = StiG^ = M^^. This is merely one (inelegant) example 
of many bulk T^iy that produce Cardassian expansion. 

We may now investigate the meaning of the values of B{n) required by Eq.(|D, where B{n) 
is the parameter in front of the new Cardassian term in Eq. (|l|) . As mentioned previously, the 
mass scale of B has units of m^~^". We find that the corresponding mass scale is very small 
for n < 1/2, is singular at n = 1/2, and then goes over to a very large value for n > 1/2. 
Specifically, for n = 2/3 and z^g = 1, we obtain B ~ 10~^^GeV~^/^ which corresponds to a 
mass scale of lO^^GeV. In the context of extra dimensions, this large mass scale turns out to 
cancel against other large numbers in such a way that it corresponds to reasonable values of 
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the energy momentum tensor in the bulk. We find that r is roughly the age of the universe 
and we have e ~ 1 for all u. Then we have 

T° ~ (10-^GeV)^ (30) 

Although this value is not motivated, it is not unreasonable. In other words, reasonable 
bulk values can generate the required parameters in Eq.(|I|). Numerical values for other 
components of T^^ are the same order of magnitude, with the exception that T04 ~ 0. For 
the case of n = 1/3, we obtain a mass scale of 10"^''^, which cancels other small numbers 
in such a way as to again require roughly Eq.(pOD to be satisfied. The form of T^u given 
in Eq.(^) is by no means unique and has been presented merely as an existence proof; we 
hope a more elegant T^^, may be found, perhaps with a motivation for the required value of 
B. 

Discussion 

We have presented Eq. (|l|) as a modification to the FRW equations in order to suggest an 
explanation of the recent acceleration of the universe. In the Cardassian model, the universe 
can be fiat and yet matter dominated. We have found that the new Cardassian term can 
dominate the expansion of the universe after z^q = 0{1) and can drive an acceleration. We 
have found that matter alone can be responsible for this behavior (but see the comments 
below). The current value of the energy density of the universe is then smaller than in the 
standard model and yet is at the critical value for a fiat geometry. Structure formation 
is unaffected before Zeq- The age of the universe is somewhat longer. The first Doppler 
peak of the cosmic background radiation is shifted only slightly and remains consistent with 
experimental results. Such a modified FRW equation may result from the existence of extra 
dimensions. Further work is required to find a simple fundamental theory responsible for 
Eq.(p. 

Questions of interpretation remain. We have said that matter alone is responsible for 
the accelerated behavior. However, if the Cardassian behavior results from integrating out 
extra dimensions, then one may ask what behavior of the radii of the extra dimensions is 
required. The Israel conditions connect the energy density on the brane to fields in the bulk. 
The required behavior of bulk fields is not transparent when one writes the modified FRW 
equation. We have found a large or small mass scale to be required, which must result from 
the extra dimensions. In principle one would like to have a complete 5-dimensional theory so 
as to perform post-Newtonian tests on the model and also to check other consequences. For 
example, with a 5-dimensional model, one would like to compare with limits from fifth force 
experiments and to check that none of the higher dimensional fields are overcontributing to 
the energy density of the universe at any point (the moduli problem). 

One might attempt to use a Cardassian term (the second term in Eq.(|I])) to drive an 
early infiationary era in the universe as well. For n < 1/2 one could have a superluminal 
expansion during the radiation dominated era. However, once accelerated expansion begins, 
Eq.(|I]) without a potential provides no way for infiation to stop. Hence we have focused 
on using this new term to generate acceleration today rather than to cause an infiationary 
mechanism early on. However, it may be possible to combine such a term with a different 
way to end inflation. 
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